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Abstract. The anisotropy of the sky distribution of 2025 
gamma-ray bursts (GRBs) collected in Current BATSE 
catalog is confirmed. It is shown that the quadrupole term 
being proportional to ^ sin 2b sin I is non-zero with a prob- 
ability 99.9%. The occurrence of this anisotropy term is 
then supported by the binomial test even with the proba- 
bility 99.97 %. It is also argued that this anisotropy cannot 
be caused exclusively by instrumental effects due to the 
non- uniform sky exposure of BATSE instrument; there 
should exist also some intrinsic anisotropy in the angu- 
lar distribution of GRBs. Separating GRBs into short and 
long subclasses, it is shown that the 251 short ones are dis- 
tributed anisotropically, but the 681 long ones seem to be 
distributed still isotropically. The 2-sample Kolmogorov- 
Smirnov test shows that they are distributed differently 
with a 98.7% probability. The character of anisotropy sug- 
gests that the cosmological origin of short GRBs further 
holds, and there is no evidence for their Galactical origin. 
The work in essence contains the key ideas and results of 
a recently published paper ( Balazs et al. 1998 ), to which 
the new result following from the 2-sample Kolmogorov- 
Smirnov test is added, too. 
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1. Mathematical considerations 

From the mathematical point of view, the necessary con- 
dition for the isotropy is the stochastic independency 
of the sky distribution of the bursts on their observed 
physical properties. This means that, if f{b,l,xi,...,Xn) 
dFdxi...dx2 is the probability of finding an object in the 
dF = cos 5 dl db infinitesimal solid angle and in the 
{xi, xi + dxi, Xn,Xn + dxn) interval, one must have 



Here Q < I < 360°, -90° < & < 90° give the celestial 
positions in Galactical coordinates, x„ (n > 1) measure 
the physical properties (peak fluxes, fluences, durations, 
etc..) of GRBs and g is their probability density. One 
may assume the fulfilment of this equation for GRBs. 

In the case of spatial isotropy, assuming that the detec- 
tion probability does not depend on the celestial direction, 
one has: uj{l,b) — l/(47r). In general case one may decom- 
pose the function uj{b, I) into the spherical harmonics. One 
obtains: 

u;{b,l) = (47r)-i/2^o- 
(3/(47r))"'^/^((jJi._i cos6sin/ — wi^i cos&cosZ + ujifl sin6)-|- 

(15/(167r))"'^^^(ijj2.-2 cos^ b sm2l + LU2.2 cos^ b cos2/— 
0^2.-1 sin 26 sin/ — 0^2,1 sin26cos/)-|- 

1) + higher order harm. (2) 



(5/(167r))i/V,o(3sin2&- 



f{l,b,xi, ...,Xn) = uj{l,b)g{xi, ...,x„). 
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The first term on the right-hand side is the monopole 
term, the following three ones are the dipole terms, the 
following five ones are the quadrupole terms. Since cu is 
constant for isotropic distribution, on the right-hand side 
any terms, except for wq, should be identically zeros. To 
test this hypothesis one has to compute the values of the 
corresponding spherical harmonic at the celestial positions 
of the observed GRBs and apply, e.g., the Student test in 
order to see that the mean of the computed values signif- 
icantly differs from zero. 

A further trivial consequence of the isotropy is the ex- 
pected equal number of bursts in celestial regions of equal 
areas. For example, one may divide the celestial sphere 
into two equal areas, e.g., taking those regions in which 
the sign of a given harmonic is either positive or negative, 
respectively. Then one may compare the number of GRBs 
in these regions by the standard binomial (Bernoulli) test. 
The details (together with the relevant references) of this 
test and also of the test based on t he spherical harmonics 
are discussed in [Balazs et al. 1998 
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Table 1. Results of the binomial test of subsamples of 
GRBs with different durations. N is the number of GRBs 
at the given subsample, kobs is the observed number GRBs 
at the first area for this subsample, and % is the proba- 
bility in percentages that the assumption of isotropy still 
holds. 



The application of the 2-sample Kolmogorov-Smirnov 
test (Press et al. 1992) on W2,-i shows that the signifi- 



sample 


N 


^obs 


(TV - k,ts) 


% 


all GRBs 


932 


430 


502 


2.0 


Tm < 2s 


251 


103 


148 


0.55 


Tgo > 2s 


681 


327 


354 


32 



2. On the existence of the intrinsic anisotropy 

In order to test the isotropy of 2025 GRBs we test the 
three dipole and five quadrupole terms. One obtains that, 
except for the terms defined by UJ2.-1 and uj2,-2, the re- 
maining six terms may still be taken to be zero. This 
means that there is a clear anisotropy defined by term 
~ sin 26 sin/. The probability that this term is zero is 
0.1%. In addition, the second quadrupole term being pro- 
portional to cos^6sin2/ is non-zero, too, with the proba- 
bility 0.6%. 

A straightforward counting of GRBs in those regions 
of equal areas, where sin2&sinZ has either positive or neg- 
ative signs, respectively, shows that 930 GRBs are in the 
first area and 1095 are in the second one. Taking p = 0.5 
probability for the binomial (Bernoulli) test, one obtains a 
0.03% probability that this detected distribution is caused 
by a chance. The observed distribution of all GRBs on sky 
is anisotropic with a certainty. 

Instrumental effects of BATSE instrument should also 
play a role, since the sky exposure of BATSE instrument 
is non-uniform. The known dependence of the detection 
probability of this instrument predicts a similar kind of 
anisotropy. Hence, the question is the following: Is this 
anisotropy caused either exclusively by the non-uniform 
sky-exposure function of BATSE instrument, or is there 
also an intrinsic anisotropy in the distribution of GRBs? 

To clarify the situation we divided the GRBs into two 
groups according to their durations Tgo. We excluded the 
dimmest GRBs. Then, 932 GRBs were separeted into the 
"short" ones (251 GRBs; Tgo < 2 s), and "long" ones (681 
GRBs; Tgo > 2 s). 

Dividing the sky into the two equal areas, as described 
above, we obtain a different behaviour for the short and 
long GRBs, respectively (see Table 1.). 

Table 1. shows that the short GRBs arc further dis- 
tributed anisotropically; there is a smaller than 1% prob- 
ability of isotropy. On the other hand, the long GRBs can 
still be distributed isotropically. Up to this point all re- 
sults mentioned in this Section, were also written down in 



cance of the difference among the samples of short and 
long GRBs is 98.7%. The short and long ones are ob- 
viously distributed differently with a probability 98.7%. 
Note here that this important quantification of the differ- 
ent behaviour of two subclasses is a new result not pre- 



sented in Balazs et al. 199i. 



We mean that these values confirm the expectation 
that there must exist some intrinsic anisotropy in the dis- 
tribution of GRBs. Once there were an exclusive instru- 
mental origin of the anisotropy all GRBs, the character of 
anisotropy should be the same for both types of GRBs; 
there should exist no difference among the short and long 
samples. Of course, the character of anisotropy is quite 
different than expected for the Galactical origin. For this 
one would need a clear non-zero 102,0 spherical harmonic. 
Hence, there is no doubt concerning the cosmological ori- 
gin of GRBs. 

3. Conclusions 

1. The distribution of 2025 GRBs is anisotropic with a 
99.97% probabihty. 

2. Separating GRBs into the two classes it is shown that 
the short ones are distributed anisotropically with a 
higher than 99% probability, but the long ones can 
still be distributed isotropically. There is a 98.7% prob- 
ability that the short and long GRBs are distributed 
differently. 

3. This different behaviour of two samples cannot be 
caused exclusively by instrumental effects; there must 
exist also some intrinsic anisotropy, too. 

4. The cosmological origin is not queried by this anisot- 
ropy. 
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